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1.  Introduction 

We  have  developed  a  bulk  PBL  model  with  a  simple  internal  vertical  structure  and  a  simple 
second-order  closure,  designed  for  use  as  a  PBL  parameterization  in  a  large-scale  model.  The 
model  allows  the  mean  fields  to  vary  with  height  within  the  PBL,  and  so  must  address  the  vertical 
profiles  of  the  turbulent  fluxes,  going  beyond  the  usual  "mixed-layer"  assumption  that  the  fluxes  of 
conservative  variables  are  linear  with  height.  This  is  accomplished  using  the  same  convective  mass 
flux  approach  that  has  also  been  used  in  cumulus  parameterizations.  The  purpose  of  this  brief 
paper  is  to  show  that  such  a  mass  flux  model  can  include,  in  a  single  framework,  the 
"compensating  subsidence"  concept,  downgradient  mixing,  and  well-mixed  layers. 


2.  Generalized  mass  flux  model 

The  "convective  mass  flux"  concept  introduced  by  Arakawa  has  been  used  in  a  variety  of 
modeling  and  observational  studies,  most  directed  either  at  deep  cumulus  convection  or  at 
stratocumulus  convection;  the  connections  between  the  two  applications  have  not  been  made 
previously. 

For  an  ascending  region  with  fracitonal  area  o,  the  conservation  equations  for  mass  and  an 
arbitrary  scalar  y  can  be  written  as 

-^(p<7)  =  -(Tv  (^pvj-Hp  (p-v)-p^(cu„ct), 
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for  the  ascending  parcels.  Similar  equations  can  be  written  for  the  descending  regions.  Here  the 
generic  variable  \\f  represents  an  intensive  thermodynamic  variable  such  as  the  potential  temperature 
or  the  the  mixing  ratio  of  water.  Other  notation  is  as  follows:  p  is  the  density;  V  is  the  horizontal 
velocity  vector,  O)  =  -pgw,  where  w  is  the  vertical  velocity,  and  g  is  the  acceleration  of  gravity;  and 
p  is  pressure.  The  air  entering  rising  parcels  has  been  assumed  to  have  the  average  properties  of 
the  sinking  parcels,  and  vice  versa.  Area  averages  satisfy 

(~)  =  (  )u^  +  i 

(2.3) 

where  an  overbar  denotes  an  area  average  (over  a  grid  box,  say);  this  notation  will  be  used  only 
where  necessary  to  avoid  confusion.  Subscripts  u  and  d  denote  upward  and  downward  moving 
parcels,  respectively.  We  can  interpret  |i.-‘  and  v*  as  time  scales  for  mass  flow  from  downdrafts 
into  updrafts,  and  from  updrafts  into  downdrafts,  respectively.  The  vertical  turbulent  flux  of 
due  to  small  eddies  is  denoted  by  fy.  The  small-eddy  fluxes  are  logically  necessary  near  the  top 
and  bottom  of  the  PBL,  where  the  organized  vertical  motions  associated  with  the  convective 
circulations  must  vanish.  The  time  scale  for  destruction  of  the  convective  circulations  by  lateral 
mixing  due  to  smaller-scale  turbulence  is  The  source  of  \\f  per  unit  mass  per  unit  time  is 
denot^  by  Sy.  This  could  represent,  for  example,  the  effects  of  rattiation. 

From  (2.1)  and  (2.2)  and  the  correponding  "downdraft"  equations,  and  using  (2.3),  we 
canrecover  the  area-averaged  conservation  equations  for  mass  and  y; 

(2.4) 

(2.5) 

Naturally,  all  terms  involving  n,  v,  and  have  dropped  out  of  (2.4-5).  In  (2.5),  the  turbulent 
fluxes  associated  with  the  convective  circulations  are  represented  by 


where 
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is  a  convective  mass  flux,  which  has  been  defined  in  such  a  way  that  it  is  less  than  or  equal  to  zero. 
The  total  turbulent  flux  of  y  is  F^  =  /,  + 


From  the  "updraft"  mass  conservation  equation  (2.1)  and  the  corresponding  "downdraft" 
equation,  we  can  derive  an  expression  for  9o/dt 
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By  neglecting  the  local  time  rate  of  change  and  large-scale  advection  terms,  we  can  simplify  (2.8) 
to  a  more  recognizable  form: 

/  dp  =  n-v. 


(2.9) 


According  to  (2.9,  the  vertical  distribution  of  (o.  is  closely  related  to  the  profiles  of  |i  and  v. 

We  can  also  derive  the  conservation  equation  for  the  variance  o/  v^that  is  associated  with 
the  convective  circulations: 
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In  (2.10),  we  have  neglected  the  local  time-rate-of-change  and  large-scale  advection  terms, 
although  we  use  them  in  our  numerical  model.  According  to  (2.12),  transport  balances  the 
combined  effects  of  dissipation  and  production.  Using  conventional  Reynolds  averaging,  we  find 
corresponding  to  (2.10)  that 


(2.12) 


Here  is  the  molecular  dissipation  rate.  Comparison  of  (2. 10)  with  (2. 12)  reveals  a  term- by- term 
correspondence.  Note  that,  in  (2, 10),  the  sign  of  the  triple  correlation  is  determined  by  the  value 
of  o.  Variance  is  transported  upwaid  for  a  <  1/2,  and  downward  for  o  >  1/2. 

The  damping  factor  a(l-o)  Xii,  which  appears  in  (2.10),  can  be  simplified  by  noting  that 
small  eddies  will  be  most  effective  at  reducing  the  differences  between  the  ascending  and 
descending  branches  of  the  plumes  when  one  of  the  two  branches  is  much  narrower  than  the  other, 
and  least  effective  when  the  two  branches  have  the  same  width.  On  this  basis,  we  assume  that  Til. 
s  0(1-0)  Tii  is  independent  of  o.  We  also  adopt  the  definition  (7^^,)-*  s  -A  0}V5pm>  where  A  is 
positive  and  nondimensional.  This  is  motivated  by  the  idea  that  the  capacity  of  small  eddies  to 
dissipate  the  convective  circulations  by  lateral  diffusion  increases  as  co.  increases,  and  as  the  PBL 
depth  decreases,  so  that  A  should  be  less  variable  than  T^i,. 


Assuming  that  the  inflow  and  outflow  layers  do  not  overlap,  we  can  eliminate  p  and  v  in 
(2.12),  to  obtain 
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The  0).  term  of  (2.14)  is  never  negative,  so  it  contributes  to  X  in  the  same  sense  as  the  positive  A 
term.  It  follows  that  X  >  0  unless  a  increases  strongly  upward.  The  O).  and  A  terms  of  (2.14) 
both  arise  from  dissipative  processes.  The  3o/9p  term  of  (2.14)  arises  from  the  transport  term  of 
(2.10). 


3.  Discussion 

Several  conclusions  can  be  drawn  from  (2.13-14),  by  considering  the  case  =  0,  so  that 

Z^  =  0. 


First,  if  (l-2a)/X  is  small  (a  near  1/2  and  strong  dissipation),  we  get 


Q). 
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This  is  a  downgradient  diffusion  formula,  with  diffusion  coefficient  -O)./  (gX). 

Next,  consider  a  layer  that  is  well-mixed  in  several  conservative  variables  \j/i,  Vz,  etc.  In 
order  to  satisfy  (2.13)  simultaneously  for  all  of  the  variables  we  must  take 


X  =  0, 

<7  =  1/2. 


(3.2) 

(3.3) 


Since,  according  to  (3.3),  a  is  a  constant,  (2.14)  implies  that  X  is  the  sum  of  two  non-negative 
terms.  Then  (3.2)  implies  that  each  term  must  be  zero,  so  that 


0).  =  constant, 

(3.4) 


A  =  0. 


(3.5) 


Referring  back  to  (2.10),  we  find  a  trivial  balance  in  which  gradient  production,  transport,  and 
dissipation  are  all  zero.  (When  a  =  1/2,  the  variance  transport  by  triple  correlations  vanishes.) 

Finally,  suppose  that  X  is  small  (weak  dissipation),  and  that  O  is  small.  Then  (2.14) 
reduces  to 


3p  ~  g  dp' 

(3.6) 

This  balance  is  characteristic  of  the  temperature  and  rooisnire  budgets  of  tropical  cumulus  layers,  in 
which  "compensating  subsidence"  produces  the  convective  effects  on  the  mean  profiles.  Arakawa 
(1969)  derived  (3.6);  it  has  also  been  discussed  by  Ooyama  ,  among  others.  [The  "detrainment" 
term  included  in  cumulus  parameterization  theories  is  also  present  in  (2.13-14);  it  arises  from  the 
gradient  of  to.  in  (2.14).] 


For  the  cumulus  regime,  the  conservation  equation  for  y  can  often  be  approximated'by  a 
balance  between  vertical  advection  and  convective  transport,  i.e. 


df^  _  CD  dy/ 
^  ~  g  dp' 


Comparing  (3.6)  and  (3.7),  we  find  that 

<T  =  (l-a>./m)/2. 


(3.7) 


(3.8) 


For  tropical  cumulus  layers  with  large-scale  rising  motion,  to.  =  o),  so  that  (3.8)  is  consistent  with 
observations  that  show  a  «  1.  For  the  case  of  large-scale  sinking  motion,  (3.8)  predicts  larger 
values  of  a,  again  in  qualitative  agreement  with  observations. 

These  results  show  that  (2.14-15)  are  consistent  with  a  variety  of  observed  balances  in 
convective  layers.  TTicy  provide  a  dynamical  interpolation  between  the  "con^nsating  subsidence" 
and  "mixing  length"  regies. 
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